Matrix elements and Parton Shower in the event generator BABAYAGA 
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A new version of the event generator BABAYAGA is presented, which is based on an original matching of the Parton 
Shower approach with the complete exact 0(a) matrix element for the inclusion of the QED radiative corrections 
to the Bhabha process at flavour factories. The theoretical accuracy of the improved generator is conservatively 
estimated to be 0.2%, by comparison with independent calculations. The generator is a useful tool for precise 
luminosity determination at flavour factories, for center of mass energies below 10 GeV. 



1. Introduction 

The precise determination of the machine lumi- 
nosity is an important ingredient for the success- 
ful achievement of the physics programme at the 
e + e~ colliders running with center of mass energy 
in the range of the low lying hadronic resonances. 

One of the most important challenges is the 
precise measure of the R ratio, by means of the 
energy scan or the radiative return method. The 
aim is to reduce the theoretical error on the 
hadronic contribution to the vacuum polariza- 
tion, which reflects on the error of the anomalous 
magnetic moment of the muon a M and the QED 
coupling constant at the Z peak cxqed{Mz) pQ. 
The precise measure of R will give a stringent test 
of the Standard Model predictions. 

In order to achieve a precise determina- 
tion of the collider luminosity, high-precision 
calculations of the QED processes e + e~ — > 
e + e _ , 77, and relative Monte Carlo gener- 

ators, are required. The large-angle Bhabha pro- 
cess is, in particular, of major interest for its large 
cross section and its clean experimental signature. 

The event generator BABAYAGA was devel- 
oped |2l:i| to this end and it is widely used by 
experimental collaborations. The event genera- 
tors MCGPJ g|, BHAGENF [5] and BHWIDE are 
also employed by experimental collaborations. 

The theoretical precision of the BABAYAGA event 
generator for the Bhabha process was estimated 



at the 0.5% level |2J. The need of an improvement 
in the calculation of radiative corrections within 
BABAYAGA emerged in the last years, for a num- 
ber of reasons. First, the total luminosity error 
quoted by KLOE is presently 0.6% 0, where the 
dominant source of uncertainty comes from the- 
ory, i.e. from the 0.5% physical precision of the 
BABAYAGA generator. Secondly, the measurement 
of the hadronic cross section in the ir + ir~ channel 
at VEPP-2M has achieved a total systematic un- 
certainty of 0.6%, requiring an assessment of the 
collider luminosity at the level of 0.1%. Last but 
not least, precision measurements of R trough ra- 
diative return at KEK-B and PEP-II are already 
performed or foreseen in the near future. 

Here, we describe the main features of a high- 
precision calculation of photonic radiative cor- 
rections to the Bhabha process, in order to im- 
prove the theoretical formulation of the original 
BABAYAGA generator up to O(0.1%). The ap- 
proach is based on the matching of exact next- 
to-leading (0(a)) order corrections with resum- 
mation through all orders of a of the leading con- 
tributions due to multiple radiation, taken into 
account according to a QED Parton Shower (PS). 

2. Matching exact 0(a) matrix element 
and Parton Shower 

The original version of BABAYAGA (3.5) is based 
on a PS in QED which allows to calculate the 
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photonic radiative corrections to the Bhabha pro- 
cess. The PS is a numerical algorithm which ex- 
actly solves the DGLAP equation for the electron 
Structure Function (SF) in QED, allowing to in- 
clude essentially the leading logarithmic (LL) cor- 
rections to the cross section up to all orders of a. 
The main advantage of the PS approach is that, 
thanks to its Monte Carlo nature, the events can 
be generated exclusively, i.e. all the momenta of 
the final state particles (fermions and an indef- 
inite number of photons) can be reconstructed. 
The PS can be improved to account also for in- 
terference effects between initial and final state 
radiation, allowing thus for a more accurate de- 
scription of the radiative events [3] . 

The corrected differential cross section in the 
PS approach can be written (in a simplified form 
for the sake of clarity) as 
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where F$v an d F^u are identically equal to 1 
in the pure PS case, U(Q 2 ,e) = exp(— ^-I + L) is 
the Sudakov form factor (Q 2 = st/u is the scale 
in the collinear logarithm L = \og(Q 2 /m 2 ) — 1 
and J+ = J P(y)dy), P(x) is the Altarelli-Parisi 
splitting function, xt is the energy fraction carried 
by the i th photon, e is an infrared separator (the 
integrated cross section is independent from its 
value) and I{ki) is the photon angular spectrum 
function including initial-final state interference. 
In Eq. the photonic radiative corrections are 
resummed up to all orders, in the LL accuracy. 
The expansion at 0(a) of Eq. is 
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where the subscript SV defines the cross section 
with the emission of a virtual or soft photon (with 
energy fraction lower than e) and H defines the 
hard bremsstrahlung cross section (with energy 
larger than e). The PS error is due to the fact 
that both the soft plus virtual and hard part are 



approximated at the LL accuracy. By comparing 
the cross section of Eq. (JSJ with the exact 0(a) 
cross section it was possible to establish that the 
theoretical accuracy of BABAYAGA was at 0.5%. 

In order to go beyond the LL approximation 
and preserve the resummation of the higher or- 
ders (h.o.) corrections, an appropriate procedure 
must be devised, able to avoid the double count- 
ing of the LL contributions and guarantee the in- 
dependence of the cross section from the e param- 
eter. A solution to this problem can be obtained 
by setting the factors F$v and Fi t H to 
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where dagy X and da^ ex are the complete ex- 
pressions for the soft plus virtual 0(a) cross sec- 
tion jS] and the real one photon emission cross 
section. By construction, the factors of Eq. © 
are infrared safe quantities and they let the 0(a) 
expansion of da°° in Eq. (JJJ coincide with the 
exact 0(a) cross section. It is worth noticing 
that correcting each single photon emission with 
a factor is crucial to make the cross section 
independent from the e parameter and that the 
correction is larger in those phase space regions 
where the PS is more unreliable, typically where 
the photon is hard and/or not collinear to one of 
the charged particles. Notice that Eq. (JIJ is cast 
in a completely differential form, so that events 
can be generated exclusively as in the pure PS 
approach. 

A more detailed discussion of the match- 
ing procedure, as well as its implementation in 
BABAYAGA, is presented in Ref. [§]. 

3. Numerical results 

In order to check the technical implementa- 
tion and the physical accuracy of Eq. (JTJ in 
the new BABAYAGA ^U], the generator has been 
compared with independent calculations, namely 
LABSPV PJ, BHWIDE and the old 3.5 release, tun- 
ing the relevant input parameters. In LABSPV, 
which is not a true event generator, the corrected 
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Table 1 

Comparison of BABAYAGA with BHWIDE, LABSPV and the old version 3.5. Cross sections are in nb. 



setup 


Born 


0(a) 


BABAYAGA 


BHWIDE 


LABSPV 


BABAYAGA 3.5 


a) 
b) 


6855.74(1) 
529.463(1) 


6059.9(1) 
451.53(1) 


6086.61(2) 
455.853(4) 


6086.3(2) 
455.73(1) 


6088.5(3) 
456.20(1) 


6107.7(2) 
458.44(1) 



cross section is calculated as a = (1 + Cnl)<^sf, 
where Cnl = (a a,ex — a'g F )/ao and the Structure 
Functions in the strictly collinear limit are used: 
thus 0(a) corrections are included exactly and 
h.o. are included in the collinear limit with LL 
accuracy. BHWIDE is an event generator based on 
the YFS formalism to exponentiate exact 0(a) 
corrections. Here, O(o) is included exactly and 
h.o. are included in the YFS approach. 

The comparisons have been performed neglect- 
ing vacuum polarization effects, at 1.02 GeV of 
center of mass energy, requiring for the final 
state e + and e~ E± > 0.4 x Ec»Mi the e + e~ 
acollinearity £ lower than 10° and choosing two 
angular acceptances: a) 20° < #± < 160° and 
b) 55° < 0± < 125°. 

The results of the comparisons on the inte- 
grated cross section are shown in Tabled where 
also the Born and 0(a) predictions are added. 
In setup a) (b)), the 0(a) corrections reduce the 
Born cross section by 13% (17%) and the cor- 
rections beyond 0(a) raise the cross section by 
0.4% (0.9%), showing that h.o. contributions are 
important at the level of the needed accuracy. 
The agreement between BABAYAGA and BHWIDE is 
very good, being the discrepancies well below the 
0.1%. The tiny differences with respect to LABSPV 
(0.08% in setup b)), where the impact of h.o. is 
larger) can be ascribed to the strictly collinear 
approximation of h.o. corrections in LABSPV. Fi- 
nally, the differences between the new and old 
version of BABAYAGA arc within the 0.5%, in agree- 
ment with the theoretical accuracy estimated for 
the old release of the event generator. 

Besides the integrated cross section, also dif- 
ferential cross sections have been compared. In 
Fig-H the acollinearity distribution is plotted for 
setup b), as obtained with the old and the new 
BABAYAGA and at 0(a). The effect of h.o. cor- 
rections is clearly visible. The relative difference 
between the old and the new BABAYAGA is better 



pointed out in the smaller panel, being almost 
constant in the whole range at 0.5% level. 
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Figure 1. Acollinearity distribution, the new 
BABAYAGA compared with its old version. 



In Fig. |3 the comparison of the acollinearity 
distribution obtained with BABAYAGA and BHWIDE 
is shown. The differences are very small and can 
be only appreciated in the panel, showing they 
reach the 0.4% level only in the distribution tail. 
This region is populated by events where hard 
photons are emitted and where the different treat- 
ment of the multiple-photon emission and h.o. 
corrections shows up, as can be expected. 

4. Estimate of the theoretical accuracy 

The results presented the previous section 
clearly demonstrate that the matching procedure 
of Eq. (JTJ has the desired behaviour of includ- 
ing the missing 0(a) contributions in the PS al- 
gorithm while preserving the resummation of all 
h.o. corrections. They also demonstrate that the 
not trivial Monte Carlo implementation of Eq. Q 
is correctly realized and robust. 
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Figure 2. Acollinearity distribution, BABAYAGA 
compared with BHWIDE. 



The theoretical error of BABAYAGA is now 
shifted at the two loop level, to terms of the or- 
der of \a 2 L (~ 3 • 1CT 4 at 1 GeV), without any 
infrared enhancement. 

In [Hj, a critical and systematic discussion of 
the remaining theoretical error will be presented. 
Once the vacuum polarization is (easily) included 
in the BABAYAGA formulation, the formulae can be 
expanded at order a 2 and consistently compared 
with the recent calculation of Refs. where 
relevant contributions to the complete two loop 
corrections to Bhabha scattering are computed. 
It is worth noticing that in Refs. the real pho- 
ton radiation is accounted only in the soft limit 
and that at two loop level also the emission of 
real e + e~ pairs has to be carefully considered. 
However, preliminary studies and tests show that 
the missing 0(a 2 ) contributions have an impact 
not larger than the 0.1% on the integrated cross 
section if typical event selection criteria for lumi- 
nometry at flavour factories are applied. 

These considerations suggest that a conserva- 
tive estimate of the theoretical error of the im- 
proved BABAYAGA can be fixed at 0.2%. 

5. Conclusions 

An improved version of the event generator 
BABAYAGA is now available ^U], where an origi- 
nal matching of the QED PS and the exact O(a) 



matrix element is implemented in order to go be- 
yond the LL approximation intrinsic to the PS 
approach. The BABAYAGA theoretical error for the 
Bhabha cross section calculation at flavour facto- 
ries is reduced from the 0.5% down to 0.2% (con- 
servatively estimated). A more detailed descrip- 
tion of the implementation and a more robust es- 
timate of the theoretical accuracy are discussed 
in Ref. 0. 
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